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Abstract. We present results of a study of the Galactic ridge X-ray emission (GRXE) in hard X-rays with the 
IBIS telescope on board INTEGRAL in the region near the Galactic Anticenter (GA) at I = 155°. We assumed 
a conservative 2a upper limit on the flux from the GA in the 25 — 60 keV energy band of 1.25 x 10 -10 erg s _1 
cm -2 (12.8 mCrab) per IBIS field of view, or 6.6 x 10 -12 erg s _1 cm -2 (0.7 mCrab) per degree longitude in the 
135° < I < 175° region. This upper limit exceeds the expected GRXE intensity in the GA direction by an order 
of magnitude, given the near-infrared (NIR) surface brightness of the Milky Way in this region and the standard 
hard X-ray-to-NIR intensity ratio for the GRXE, assuming stellar origin. Based on the CGRO/EGRET surface 
brightness of the Galaxy above 100 MeV as a tracer of the cosmic-ray (CR) induced gamma-ray background, the 
expected GRXE flux in GA exceeds the measured 2a upper limit by a factor of 4. Therefore, the non-detection 
of hard X-ray emission from the GA does not contradict the stellar nature of the GRXE, but is inconsistent with 
CR origin. 
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1. Introduction 

The stellar origin of the Galactic hard X-ray back- 
ground, better known as the Galactic ridge X-ray 
emission (GRXE), has recently been strongly sup- 
ported by morphological/spectral studies with the RXTE 
and INTEGRAL observatories (iRevnivtsev et al~ l2006t 
iKrivonos et al. I l2007at iTiirler et all 1201017 spectral stud- 
ies with 
200: 



Suzaku 
and direct 



Ebisawa et al. 



X-ra y 



source 
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counts 



Yamauchi et al 
with Chandra 



(|Revnivtsev et al.ll2009l l201lh . The GRXE does not arise 
from the interaction of cosmic rays with the interstellar 
medium, as was believed before, but is associated with 
the (predominantly old) stellar population of the Galaxy, 
namely with hard X-ray emission from accreting white 
dwarfs and coronary active stars. It was dem onstrated 



(Rev nivtsev e t al. 2006; K rivonos et al. l2007ah that the 



GRXE intensity closely follows near-infrared (NIR) sur- 



face brightness over the Milky Way, which is a known 
tracer of stellar mass. 



Galactic diffuse emission in a regi on far away from the 
Galactic Center (GC) was studied bv lWorrall et al.l (|1982l ) 
with the A2 experiment on board the HEAOl satellite. 
The observed 2—10 keV emission was consistent with an 
radial exponential disk with a half-thickness of ~ 240 pc. 
It was pointed out that unresolved emission likely comes 
from discrete point sources and does not have a diffuse 
origin. In the hard X-ray domain GRXE was studied by 
Skibo et al.l (|l997l ) using OSSE observations at longitude 
I = 95°. The observed emission between 50 and 600 keV 
was suspected to contain a significant contribution from 
bright discrete sources because of the wide collimated field 
of view (~ 11°. 4 x 3°. 8), but a major part of the de- 
tected flux was interpreted as inters tellar emission from 

20001: 



* Based on observations with INTEGRAL, an ESA project 
with instruments and science data centre funded by ESA 
member states (especially the PI countries: Denmark, France, 
Germany, Italy, Switzerland, Spain), Czech Republic and 
Poland, and with the participation of Russia and the USA 



non-thermal electrons (se e also e.g. IValinia et al 



Valinia fc Tatischeff 1200 ll) . Today, thanks to the unique 
possibilities of INTEGRAL gamma-ray telescopes, we can 
directly study the Galactic diffuse background in any parts 
of the Galaxy without dealing with significant source con- 
tamination. 
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Given the distribution of 4.9/im NIR intensity over the 
Galaxy as measured by the COBE/DIRBE experiment, 
the GRXE is not expected to be detectable in the Galactic 
Anticenter (GA) because of the low NIR surface brightness 
in this region. Nevertheless, an explicit demonstration that 
the GRXE is not observed from a Galactic region of low 
stellar density, such as the GA, would substantiate the 
stellar paradigm of the GRXE even more. Placing tight 
constraints on the hard X-ray flux from the GA region is 
also important for calibrating future studies of the GRXE 
in the central parts of the Galaxy. 

2. Observations 

We used data from the ISGRI detector, the first 



Table 1. INTEGRAL observations used for the GRXE 
study. 



layer of the IBIS coded-mask telescope (jUbertini et al 



2003) , on board th e INT EGRAL gamma-ray observa- 
tory (jWinkler et al.l 120031 ). ISGRI operates at energies 
above ~ 20 keV, with the sensitivity rapidly decreasing 
above 100 keV. IBIS has a relatively wide field of view 
(~ 28° x 28° at zero response), which allows one to mea- 
sure weak diffuse emission fluxes by using the telescope as 
a collimated instrument. 

To study the GRXE in the GA region, we used spe- 
cial INTEGRAL observations, part of a series of so-called 
Galactic latitude scans (GLS). This program is based on 
consecutive observations made along the Galactic latitude 
in the range ±30° with a step of 2°. This strategy al- 
lows one to make independent snapshot measurements of 
the instrumental background at mid latitudes, where the 
GRXE is expected to be negligible, along with an actual 
GRXE observation near the Galactic plane. It is crucial 
that the instrumental background, which is usually high 
in hard X-ray observations, does not exhibit strong vari- 
ability during an individual GLS lasting ~ 8 hours. 

There are two positions in the GA region observed with 
GLSs: ongoing observations at I — 215° (PI: Tsygankov) 
and a completed program at I — 155° with a total ex- 
posure of 1 Ms (PI: Krivonos). In the present study we 
only used the completed GLS observations at I = 155° 
performed in August-September 2010, see Table [T] for de- 
tails. After screeni ng the whole data set fo llowing the pro- 
cedure described in Krivonos et al. I (<2007ah , hereafter K07, 
we were left with 525 out of 569 short (~ 2 ks) individ- 
ual INTEGRAL observations, called scientific windows 
(ScWs), for the subsequent scientific analysis. To model 
background variations we used public data of high-latitude 
observations (Table [lj. 

3. Analysis 

We mainly followed the approach described in K07. Our 
study of the GRXE is based on the capability of the IBIS 
telescope to separate the contributions of point sources 
from the background count rate. 

The ISGRI detector shadowgram that was accumu- 
lated during individual INTEGRAL observation in a given 
energy range was cleaned from source fluxes using the 



Latitude scan at I = 155°, 2010 


Observation field Orbits 


ScWs 


960,961 


2-94,11-97 


962,963 


12-112,1-97 


964,965 


13-110,2-67 


966 


23-56 


569 ScWs, total nominal exposure: 1 Ms. 


Background model calibration, 


2008-2009 


Observation field Orbits 



North Ecliptic Pole 
XMM-LSS 
Virgo Cluster 
Coma Cluster 
M82 X-l 



686-689, 759-761, 824-829 

695, 696, 701 

747-754, 758, 819-820 

821-823 

869-872 



3C273, 3C279 and M87 878-880 



Total nominal exposure: 6 Ms. 



Crab calibration, 2010 



Observation field 



Orbits 



Crab Nebula 



902, 903, 966-968, 970 



Total nominal exposure: 890 ks. 



known pattern of the IBIS mask. The remaining shad- 
owgram contains the following components: (i) variable 
instrumental background, (ii) the isotropic cosmic X-ray 
background (CXB), and (iii) the GRXE. Throughout the 
analysis we assumed the CXB flux as a constant part of 
the instrumental background. The possible influence of the 
CXB variance on the GRXE measurement is considered in 
Appendix [B] Because the GRXE is extended over the sky, 
it cannot be directly resolved with the IBIS mask. The 
only way to estimate the GRXE flux is to determine the 
difference between the observed collimated detector count 
rate that is cleaned from sources and the assumed instru- 
mental background. To this end, we define a background 
model in Appendix IA1 that predicts the background count 
rate during a given GRXE observation. The background 
model is adjusted very precisely using mid-latitude snap- 
shots of the background performed shortly before and af- 
ter a given GRXE observation, which is the main concept 
of the GLS. 



3.1. Detector shadowgram 

Because IBIS is a coded aperture imaging telescope, the 
sky is projected onto the detector plane through the trans- 
parent and opaque elements of the mask mounted above 
the detector plane. We produced the ISGRI detector shad- 
owgram for every ScW as described in K07. We used the 
25 — 60 keV working energy band because of the known 
evolution of the low-energy threshold of the ISGRI detec- 
tor and because the GRXE is expected to be weak above 
~ 50 keV owing to the high-ene r gy cut-off in the GRXE 
spectrum ( Krivonos et al. 2007at Revnivtsev et al. I l2006lb 
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3.2. Sky map 

The sky reconstruction is based on deconvolution of the 
detector shadowgram with a known mask pattern. We 
implemented the IBIS/ISGRI sky recon struction method 
descr i bed in our previous publications (IRevnivtsev et al, 



2004: Krivono s et all 1200 



pubncatioi 
I l2007alibh 



we ref er th e reader to the paper s by 



For the basic idea 



Fenimore & Cannon 



(1981) and Skinner et al. ( 1987T ) . Every sky image was ad- 
ditionally cleaned from systematic noise as described in 
Krivonos et al.l ( 2010a ). The resulting sky image mosaic 



is shown in Fig. [TJ 

The survey area as a function of flux for sources with 
S/N > 5 is shown in Fig. [5J The minimum detectable 
flux in the central part of the field is 2.4 x 10~ n erg 
s _1 cm~ 2 (or 2.5 mCrab). The survey area reaches its ge- 
ometric limit of 3650 deg 2 
cm 



for / > 4.9 x 10" 9 erg s" 



2 (500 mCrab), 50% of this area has a sensitivity bet- 

.3 mCrab). 



ter than 8.1 x 10 n erg s 1 cm 2 



We performed a search for sources as excesses in the 
sky mosaic (Fig. [IJ convolved with a Gaussian represent- 
ing the effective instrumental PSF. The detection thresh- 
old was estimated assuming Gaussian noise of the pixel 
values as follows. The total area of the sky image (Fig.[T|) is 
3650 squared degrees, and taking the IBIS telescope angu- 
lar resolution of 12' into account, we gathered ~ 9.12 x 10 4 
independent pixels. With this consideration, we set the 
source detection threshold to 5<j, allowing at most one 
false detection for pure Gaussian noise. 

The signal-to-noise ratio distribution of pixels shown 
in Fig. [3] has the expected Gaussian shape. However, 
one notices some systematic excess at the negative side. 
Nevertheless, the source detection threshold of 5er, esti- 
mated above, separates the noise and source dominated 
pixel domains. 

The list of point sources is presented in Table [5J We 
attribute two marginally detected, previously unknown 
sources to the systematic noise. IGR J08135+5655 is 
not detected in the sky mosaic of the slightly broader 
energy band 20 — 60 keV, and the region around 
IGR J01532+2612 is affected by systematic noise at the 
edge of the sky mosaic. 

The list of sources in Table [2j except for two IGR's 
mentioned above, was used in the iterative source re- 
moval (IROS) procedure th at we applied to every ISGRI 



detector shadowgram (e.g. iKrivonos et al.l l2010a[) . This 



procedure introduces additional uncertainty to the back- 
ground model (Appendix [C]), but allows one to trace 
source variability. For example, the known Be/X-ray bi- 
nary RX J0440.9+4431 (LS V +4 4 17) was in a strong 



outburst during the o bservations ([Krivonos et al.l 12010c ; 



Tsveankov et al.l 120111 ). Fig. g] shows 25 - 60 keV detec- 
tor light curve, cleaned from the source contribution, and 
ready for further analysis. 
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Fig. 1. IBIS/ISGRI sky image of the 1 Ms observation of 
the GA region at I — 155° produced in the 25 — 60 keV 
energy range. The color table of the image represents pixel 
values in the range 0-25 calculated as the square root of 
the significance. The cataloged and newly detected sources 
are labeled in yellow and green, respectively. 



4. Results 

Using the background model (Appendix [X]) we obtained 
for each Sc W the predicted detector count rate, as shown 
in Fig. [5] The detailed view of the observed and predicted 
detector light curves and their residuals during a space- 
craft orbit 964 is shown in Fig. [5] We denoted regions of 
background measurement (|6| >= 20°) in blue, and actual 
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Fig. 4. Detector count rate of the individual 5c Ws after subtracting flux contribution from point sources during the 
GLS program at I — 155°. The data points contain pure statistical errors. The labeled horizontal intervals denote 
different spacecraft orbits. 
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GRXE observations < 20°) in red. Fig. M shows that 
the background behavior in an individual orbit can be cap- 
tured only with fast scanning observations such as GLSs. 
The background model, with low intrinsic scatter, exactly 
follows the observed detector rate, and the scatter of resid- 
uals (lower panel of Fig. [B]) is comparable to the statistical 
uncertainty of the IROS procedure (Appendix [C]) . 

Using the entire GLS data set at / = 155°, we averaged 
residuals over Galactic latitude, as shown in Fig. [7] by red 
points. The latitude profile does not show any significant 
excess in the Galactic plane region at b = 0° . As expected, 
the GRXE associated with the old stellar population is 
not detected in the GA. The la upper limit on the GRXE 
flux in the |6| < 5° latitude range, which roughly corre- 
sponds to the IBIS fully coded FOV, is ~ 2.8 mCrab, or 
~ 6.4 mCrab taking systematic uncertainties into account. 
We later refer to a 2a upper limit of ~ 12.8 mCrab. One 
can convert the achieved upper limit to more convenient 
units using the effective solid angle of the IBIS telescope 
~ 286 deg 2 and taking into account that the GRXE is 
much less extended in the Galactic latitude than the aver- 
age cross-section of IBIS FOV. For instance, the 2a upper 



.0e + 02 ¥ 




Signal — to — noise ratio 

Fig. 3. Signal-to- noise ratio distribution of pixels in the 
1 Ms GLS observation at I = 155°. The dashed line repre- 
sents the normal distribution with unit variance and zero 
mean. The accepted threshold of source detection (5<r) is 
shown by the dotted line. The plot is truncated at a = 9 
for convenience. 

limit on the GRXE flux at / = 155° per unit Galactic 
longitude is 0.7 mCrab deg -1 . 

One can test the GRXE non-detection in the GA for 
consistency with stellar and truly diffuse GRXE origins. 
To this end, we compared the observed drop of the hard 
X-ray flux from the GC to the GA region (per IBIS FOV) 
with the corresponding change of the intensity of a given 
tracer. We used the COBE/DIRBE 4.9/im datfQ as a 
tracer of stellar mass, and the EGRET gamma-ray back- 
ground map above 100 MeV as a tracer of the cosmic- 
ray induced gamma-ray background. The EGRET back- 
ground intensity drops from the GC to the GA by a fac- 
tor of ~ 3. Therefore, using the conservative estimate of 
the GRXE flux in the GC from K07 of 150 ± 15 mCrab, 

1 COBE/DIRBE 4.9/im intensity map was corrected for the 
interstellar reddening and mean background level measured in 
high-latitude regions, as described in K07. 
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Fig. 5. Detector count rate (black) of the individual ScWs as shown in Fig. [4] The systematic uncertainty of the 
source removal procedure (Appendix [Cj was added to the statistical errors of each point. The red points represent the 
count rate predicted by the background model. 



Table 2. The list of sources significantly (> 5c) detected 
on the sky mosaic (Fig. [T|). The newly detected sources 
are highlighted in bold. The 68% confidence interval for 
estimated sky coordinates depends on source significance: 
2.1', 1.5', and < 0.8'for 5 - 6,10, and > 20a, respectively 
((Krivonos et al.ll2007bh 



Name 


Ion. 


lat. 


^25-60 kcV 




deg. 


deg. 


mCrab 


4U 0352+30 a 


163.08 


-17.14 


39.42 ± 0.69 


RX J0440.9+4431 a 


159.82 


-1.26 


26.12 ± 0.53 


3Clll b 


161.68 


-8.83 


5.32 ± 0.60 


MCG 8-ll-ll b 


165.74 


10.41 


7.27 ± 0.86 


4C 32.14 c 


158.99 


-18.77 


3.15 ± 0.56 


LEDA 168563 b 


157.26 


3.43 


2.76 ± 0.49 


MRK 3 b 


143.29 


22.72 


7.23 ± 1.33 


Perseus CL d 


150.58 


-13.25 


3.06 ± 0.58 


IGR J01532+2612 


139.82 


-34.64 


20.31 ± 3.91 


IGR J08135+5655 


160.75 


33.53 


7.09 ± 1.38 


GK Per e 


150.96 


-10.12 


2.82 ± 0.56 



HMXB, b AGN, c QSO, d Cluster of Galaxies, e CV 



we expect the corresponding flux from the GA to be 
~ 50 mCrab. This is definitely not observed according 
to Fig. [7] In contrast, there is a factor of 270 drop in the 
NIR 4.9/xm intensity from 2.7 x 10~ 5 erg s^ 1 citT 2 to 10~ 7 
erg s _1 cm~ 2 . This implies a GRXE flux from the GA of 
0.4 mCrab at 25 — 60 keV, which is consistent with the 
derived upper limit of 12.8 mCrab. This is illustrated in 
Fig. El where the COBE/DIRBE and EGRET longitude 
profiles are renormalized to the hard X-ray flux observed 
from the GC. We conclude that the non-detection of the 
GRXE from the GA is consistent with the stellar mass 
distribution of the Galaxy traced by NIR maps, rather 
than with the cosmic-ray induced gamma-ray background 
seen by EGRET. 
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Fig. 6. Upper panel: Detailed view of detector count rate 
(black) during orbit 964 from Fig. [5j The blue and red ar- 
eas denote observations made at |6| >= 20° and |6| < 20°, 
respectively. Bottom panel: Residuals after subtracting the 
model-predicted count rate from the observed count rate. 
The black dashed lines represent a la deviation (1.02%) 
from zero. 



5. Conclusions 

1) Using the 1 Ms observations of the GA at I — 155° 
with the INTEGRAL observatory, performed in the spe- 
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Galactic latitude, deg. 

Fig. 7. Source flux contribution (blue) and ISGRI detec- 
tor background count rate residuals (red) averaged over 
the Galactic latitude. Error bars of the blue histogram rep- 
resent rms-deviations of the summed point source fluxes 
from average in bin. 




Galactic longitute, deg. 



Fig. 8. Galactic longitude profiles of the COBE/DIRBE 
4.9/im intensity (dashed line) and EGRET background 
above fOO MeV (dotted line), both normalized to the 
GRXE flux of 150 mCrab (red point) in the GC. The 
2cr upper limit corresponds to the GRXE measurement at 
I = 155° in the present study. 



cial GLS mode, we did not detect the GRXE in the 
25 — 60 keV energy band and set a conservative 2a upper 
limit of 1.25 x 10~ 10 erg s" 1 cm- 2 (12.8 mCrab) per IBIS 
FOV or 6.6 x 10~ 12 erg s" 1 cm" 2 deg~ x (0.7 mCrab deg" 1 ) 
per unit Galactic longitude. 

2) The obtained upper limit is consistent with the con- 
siderable drop in the NIR (4.9/i) intensity observed by 
COBE/DIRBE and disagrees with the much smaller de- 
crease in the gamma-ray (above 100 MeV) background 
measured by EGRET. Therefore, the non-detection of the 
GRXE in the GA is consistent with the stellar mass distri- 
bution in the Galaxy, which does not contradict the stellar 
nature of GRXE, but is inconsistent with its CR origin. 



3) The developed background model potentially allows 
one to reach the statistically limited accuracy. However, 
the final uncertainty of the approach is associated with the 
source removal procedure, the systematic uncertainty of 
the method itself, and the CXB variance. Nevertheless, the 
implemented method along with the special GLS mode of 
observation is an optimal approach of modeling the ISGRI 
background and can be efficiently used for studying the 
Galactic hard X-ray background. 

Acknowledgements. This research was made possible thanks 
to the unique capabilities of the INTEGRAL observatory. 
The data used were obtained from the European and Russian 
INTEGRAL Science Data Centers. The work was supported 
by the President of the Russian Federation (through the pro- 
gram of support of leading scientific schools, project NSH- 
5069.2010.2, and grant MD-1832.2011.2), by the Presidium of 
the Russian Academy of Sciences/RAS (the program "Origin, 
Structure, and Evolution of Objects of the Universe"), by 
the Division of Physical Sciences of the RAS (the program 
"Extended objects in the Universe", OFN-16), by the Russian 
Basic Research Foundation (grant 10-02-00492- A), State con- 
tract 14.740.11.0611, and the Academy of Finland grant 
127512. 

References 

Cfmrazov, E., et al. 2007, A&A, 467, 529 
Ebisawa, K., et al. 2008, PASJ, 60, 223 
Fabian, A. C, & Barcons, X. 1992, ARA&A, 30, 429 
Fenimore, E. E., Cannon T. M., 1981 Applied Optics, 20, 
1858. 

Skinner, G. K., Ponman, T. J., Hammersley, A. P., & 

Eyles, C. J. 1987, Astroph.Sp.Sci., 136, 337 
Gruber, D. E., Mattcson, J. L., Peterson, L. E., & Jung, 

G. V. 1999, ApJ, 520, 124 
Hajdas W., Biihler P., Eggel C, Favre P., Mchedlishvili 

A., Zehnder A., 2003, A&A, 411, L43 
Krivonos, R., Vikhlinin, A., Churazov, E., Lutovinov, A., 

Molkov, S., & Sunyaev, R. 2005, ApJ, 625, 89 
Krivonos, R., Revnivtsev, M., Churazov, E., Sazonov, S., 

Grebenev, S., & Sunyaev, R. 2007, A&A, 463, 957 
Krivonos, R., Revnivtsev, M., Lutovinov, A., Sazonov, S., 

Churazov, E., & Sunyaev, R. 2007, A&A, 475, 775 
Krivonos, R., Revnivtsev, M., Tsygankov, S., Sazonov, S., 

Vikhlinin, A., Pavlinsky, M., Churazov, E., & Sunyaev, 

R. 2010a, A&A, 519, A107 
Krivonos, R., Tsygankov, S., Revnivtsev, M., Grebenev, 

S., Churazov, E., & Sunyaev, R. 2010b, A&A, 523, A61 
Krivonos, R., Tsygankov, S., Lutovinov, A., Turler, M., & 

Bozzo, E. 2010c, The Astronomer's Telegram, 2828, 1 
Revnivtsev, M., Sunyaev, R., Varshalovich, D., et al. 2004, 

Astronomy Letters, 30, 382 
Revnivtsev, M., Sazonov, S., Gilfanov, M., Churazov, E., 

& Sunyaev, R. 2006, A&A, 452, 169 
Revnivtsev, M., Molkov, S., & Sazonov, S. 2008, A&A, 

483, 425 

Revnivtsev, M., Sazonov, S., Churazov, E., Forman, W., 
Vikhlinin, A., & Sunyaev, R. 2009, Nature, 458, 1142 



Krivonos et al.: Constraints on the ridge emission from the Milky Way anticenter 



7 



Revnivtsev, M., Sazonov, S., Forman, W., Churazov, E., 

& Sunyaev, R. 2011, MNRAS, 414, 495 
Skibo, J. G., et al. 1997, ApJ, 483, L95 
Tsygankov, S., Krivonos, R., Lutovinov, A., 2011, 

MNRAS, submitted 
Tiirler, M., Chernyakova, M., Courvoisier, T. J.-L., 

Lubihski, P., Neronov, A., Produit, N., & Walter, R. 

2010, A&A, 512, A49 
Ubertini P., et al., 2003, A&A, 411, L131 
Valinia, A., Kinzer, R. L., & Marshall, F. E. 2000, ApJ, 

534, 277 

Valinia, A., & Tatischeff, V. 2001, Exploring the Gamma- 
Ray Universe, 459, 153 

Winkler C, et al., 2003, A&A, 411, LI 

Worrall, D. M., Marshall, F. E., Boldt, E. A., & Swank, 
J. H. 1982, ApJ, 255, 111 

Yamauchi, S., Ebisawa, K., Tanaka, Y., Koyama, K., 
Matsumoto, H., Yamasaki, N. Y., Takahashi, H., & 
Ezoe, Y. 2008. larXiv:0810.0317l 



Appendix A: Background model 

The background model used in K07 was slightly changed 
in this study. Since the GLS observations were performed 
over a relatively short time period, we removed the long- 
term time part from the equation. Instead of using the 
gain parameter to trace orbital modulations of the back- 
ground rate, we used the spacecraft orbital phase P in 
a quadratic polynomial form. Hereafter, we consider only 
the detector count rates after removal of the contribution 
of point sources. The model of the detector background 
25 — 60 keV count rate, D^gd, is made of a linear combi- 
nation of the 600 — 1000 keV detector count rate, H, and 
phase 



Dbgd 

coefficient 



const + aH + b x P + b 2 P 



(A.l) 



The coefficient a was calculated using observations 
pointed away from the Galactic plane > 20°) where 
the GRXE is not expected to be observed (see Table [T]). 
The constant term and b coefficients were determined in- 
dividually for each spacecraft orbit from the observations 
at |6| > 20°, thus adjusting the model to the current 
background level. In fact, the constant term in Eq. IA.1I 
contains contribution from CXB and unknown intrinsic 
detector background. The last is also variable, which is 
traced by variability of this constant with time. In the 
current observations, its absolute value varies in the range 
of 800 — 900 mCrab from orbit to orbit, while permanent 
CXB contribution is expected to be at level of 550 mCrab 
(Appendix [B]). Finding the difference between the ob- 
served and predicted by Ea. lA.ll detector count rate should 
yield the GRXE excess in the Galactic plane. 

The detector count rate was converted to the conve- 
nient units of Crab flux, with 1 mCrab = 7.25 x 10~ 6 cts/s 
in the 25 - 60 keV band per IBIS FOV. The conver- 
sion coefficient was determined from observations of the 
Crab Nebula in 2010 (Table []}. A flux of 1 mCrab in 



the 25 — 60 keV energy band corresponds to 9.7 x 10 12 
erg s _1 cm~ 2 for a source with a Crab-like spectrum, 
10.0 x Edt, 1 phot cm 



Appendix B: CXB cosmic variance 

The CXB emission coming from the population of unre- 
solved extragalactic sources (active galactic nuclei, AGNs) 
is subject to Poissonian variations in the number of 
sources, intrinsic s ource variability, and nea rby large-scale 
structure (see e.g. iFabian fc Barcons 1992tl . Here, we es- 
timate the systematic limitations to the measured GRXE 
flux caused by CXB variations. 

Using the e xtraga lactic log N-log S relation from 
\ Crivonos et alJ ( 2010bl ) and following iRevnivtsev et all 
( 20081 ) , we estimated the relative uncertainty of the CXB 
flux in the 25 — 60 keV band as 



(61, 



CXB 



CXB 



5.5 x 10 2 5 , n / aXi 



n~ 1/2 

11 s 'dog J 



(B.l) 



where S maXl ii is the maximum flux of undetected sources 
in units of 10 -11 erg s _1 cm -2 , and fldeg ~ 286 is the 
effective solid angle of the IBIS telescope. We adopted the 
CXB intensity to be equal to 1.89 x 10 -11 erg s" 1 cm" 2 



deg , based on the CXB spectrum model of lGruber et al 



([19991 ) and the ~ 10% higher norma lization measured 
by INTEGRAL (|Churazov et al.ll2007l ). Using the limit- 
ing flux of the survey Eq. IB. II yields a CXB variance at 
the level of ~ 0.4%. The absolute value is 2.3 mCrab as- 
suming a - 550 mCrab CXB flux per IBIS FOV. In the 
current work, we consider the CXB variance for an area 
of three IBIS FOVs, which approximately corresponds to 
the effective area of the GA survey. 

Appendix C: Uncertainty in GRXE measurements 

The uncertainty of the background model, i.e. the accu- 
racy of the ISGRI background rate prediction, is subject to 
statistical and systematical errors. The former can be eas- 
ily estimated from the total number of counts (~3x 10 5 ) 
in the 25 — 60 keV energy band per typical ScW (~ 2 ks). 

The additional statistical effect is related to the IROS 
procedure, when the count rate attributed to a given 
source is removed from the detector using the known aper- 
ture function of the mask. To a first approximation, the 
total number of counts, S, associated with a given source 
is determined as the difference between the number of 
counts, D\, in the detector pixels illuminated by the source 
through the mask and the number of counts, Dq, in the 
detector pixels blocked by the mask: S = D\ — Dq, The to- 
tal flux on the detector is D — Dq + D\. Thus, subtracting 
the contribution of the source yields the residual detector 
flux D' = D — S = 2Dq- For a weak source in the center 
of the field of view, Dq ps D\ rj 1/2 X D and therefore the 
relative statistical uncertainty of measuring the detector 
count rate &=■ = \/2-R= increases by a factor of y/2. In 

y 2JJq V-D 

practice, a more complicated model of a point source (see 
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Fig. C.l. Relative root- mean-squared value of the detec- 
tor rate as a function of number of sources in the FOV. 



K07), implemented in the source subtraction algorithm, 
causes an even larger increase of statistical uncertainties. 

To estimate the IROS induced uncertainty as a func- 
tion of number of sources in the FOV, we studied a set 
of 22 consecutive observations without cataloged and de- 
tected sources (orbit 973, ScWs 58 — 80). The relative 
standard deviation of detector count rate, as a function of 
the number of simulated sources, N src , is shown in Fig lC.ll 
The first point at N src = 0, RMS = 0.3%, reflects the rel- 
ative statistical scatter of the data. As seen from Fig IC.ll 
the scatter rapidly increases with inclusion of sources in 
the FOV. A typical scatter of ~ 1.0% on the ScW time 
scale is observed in real data, as demonstrated in Sect.|U 

To estimate the systematic uncertainty of the method, 
we defined coefficients of the background model in Eq. IA.1I 
using the high-latitude observations (see Table [T]) in the 
South Galactic hemisphere (700 ks) and applied it to the 
North (5.3 Ms). The non-existent GRXE flux was aver- 
aged over a given INTEGRAL orbit divided into the three 
equal intervals having three parts: one in the middle and 
two adjacent, see Fig. IC.2l for reference. The middle part 
of each interval (in red) was supposed to have GRXE flux, 
and the neighboring parts (in blue) were used to correct 
the constant term (Ea lA.ip . This set-up mimics the GLS 
pattern of observations. 

The standard deviation of residuals from zero repre- 
sents the systematic uncertainty of our background model, 
which is found to be ~ 4.0 mCrab. The 25 — 60 keV de- 
tector count rate increased from ~ 2.2 to ~ 2.5 Crab over 
the considered time period, hence the relative accuracy of 
the model is ~ 0.17% of the observed background rate. 

We summarize all the discussed uncertainties related 
to the GRXE measurements in the 25 — 60 keV energy 
band. The values below are presented with respect to 
the background rate, which is assumed to be 2.5 Crab. 
Statistical uncertainties: 



Fig. C.2. Residuals after subtracting the model-predicted 
count rate from the observed count rate for the 5.3Ms 
high-latitude observations averaged over spacecraft orbital 
phase. The black dashed lines represent a la deviation 
(4.0 mCrab) of the averaged values from zero. Blue and 
red regions denote different phase intervals used for back- 
ground model calibration and actual measurements, re- 
spectively. 

— 0.30% (7.5 mCrab) - count statistics, observed in a 
typical ScW without any sources in the FOV. Since 
this value differs from the expected, it cannot be 
fully statistical. Some unexplored systematics or back- 
ground variability can contribute to the scatter of the 
observed detector count rate. 

— 1.00% (25 mCrab) - observed in a typical ScW with 
several sources in the FOV, related to the IROS pro- 
cedure. 

These uncertainties constitute the error of measurement 
and, naturally, decrease with increasing exposure. For ex- 
ample, the largest error of 25 mCrab decreases to 1 mCrab 
for a total exposure of 1 Ms. Systematic uncertainties: 

— 0.16% (4.0 mCrab) - root-mean-squared residuals af- 
ter background model subtraction from the observed 
count rate (Fig. l€T2"|) . 

— 0.09% (2.3 mCrab) - CXB variance per IBIS FOV, 

— 0.16% (4.0 mCrab) - CXB variance per GA survey 
area. 

The systematics quadratically add to the error of the mea- 
surement. The total systematic error owing to the back- 
ground model and CXB variance (per GA survey area) is 
5.7 mCrab. 



— 0.18% (4.5 mCrab) - count statistics, expected for 3 x 
10 5 counts per typical ScW (~ 2 ks), 



